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Abstract 


We study the dynamics of four-qubit W state under various noisy environments by solving 
analytically the master equation in the Lindblad form in which the Lindblad operators correspond 
to the Pauli matrices and describe the decoherence of states. Also, we investigate the dynamics of 
the entanglement using the lower bound to the concurrence. It is found that while the entanglement 
decreases monotonically for Pauli-Z noise, it decays suddenly for other three noises. Moreover, by 
studying the time evolution of entanglement of various maximally entangled four-qubit states, we 
indicate that the four-qubit W state is more robust under same-axis Pauli channels. Furthermore, 
three-qubit W state preserves more entanglement with respect to the four-qubit W state, except 
for the Pauli-Z noise. 
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I. INTRODUCTION 


Q. 


The theory of open systems has an important role in quantum information science 
The interaction between the system and the environment which leads to the relaxation and 
decoherence is the origin of the classical states. Therefore, studying open quantum systems 
can help us to protect the quantum states using proper methods [2-7], The application 
of quantum dynamics of open systems ranges from quantum cosmology to quantum op¬ 
tics and to quantum information js 15]. For instance, the detailed coherent control over 
components of quantum systems is one of the greatest challenges in quantum information 
processing. Among these efforts, the control of entanglement against the dissipative effects 
of the environment is an important issue. 

The dynamics of open quantum systems is usually described by the quantum Markov pro¬ 
cess as a semigroup of completely positive dynamical maps and the corresponding quantum 


master equation in the Lindblad form 


16 


17]. However, there are many quantum systems 


which display non-Markovian behavior in which correlations can give rise to memory ef¬ 


fects due to the environment 


18^20). The Lindblad equation as the most general Markov 


evolution equation for a density operator is given by [16] 

% = -{[Hs,P] + \(\L,pL') + {Lp,L']) , 


(1) 


where the Lindblad operator L represents the influence of the environment. 

Since the entanglement is known as a resource for performance enhancements in quantum 
information processing, understanding and control of the entanglement in open systems have 
many applications in quantum computation and quantum information. As the first attempt 
to study the dynamics of entanglement in open systems, the evolution of entanglement has 
been studied using the state changes rather than computing the entanglement of an arbitrary 


state by an entanglement measure 


2lj. The first real entanglement dynamics in open systems 


investigated the dynamics of entanglement for a pair initially entangled harmonic oscillators 
under the action of local environments {zj. Also, the evolution of the average entanglement 
of formation of random bipartite states has been analyzed in Ref. 23]. In these works, 
the authors showed that the entanglement vanishes at finite times and coherence decayed 
asymptotically at long times. 

The first study of dynamics of multipartite systems explored multiqubit systems under 
the influence of single-qubit depolarization and examined the robustness of three- and four- 
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qubit GHZ, W, and Dicke states and found that GHZ states are more robust than other 


generic states 


24], In Ref. 25], in addition to depolarizing noise, the effects of Pauli noises 


are studied on three-qubit GHZ and W states and discovered that three-qubit GHZ state 
preserves more entanglement than three-qubit W state. The dynamical evolution of N- 
qubit GHZ and W states for 2 < N < 7 coupling to the independent dephasing and thermal 
baths at zero temperature and infinite temperature has been compared numerically using 
multipartite concurrence as a quantifier of the entanglement. It is observed that the decay 
rate increases with N for the GHZ state. For the W state this phenomenon happens only 
for infinite temperature environment. Also, it is size-independent for dephasing and zero- 
temperature thermal reservoirs 26]. In Ref. 27], the robustness of 77-qubit GHZ state as 
resources for teleportation against Pauli channels and dephasing channel has been studied. 
It is found that three-qubit GHZ state is more robust than 77-qubit (N > 3) GHZ states 
under most noisy channels. It is also shown that three-qubit W state is more robust than 
three-qubit GHZ state for small noi sy p arameter while the GHZ state becomes more robust 
when the noisy parameter is large 28]. Moreover, for W state-like superpositions against 
dephasing and amplitude damping channels, using three measures of entanglement, it is 
shown that the effects of decoherence on the fidelity is better described by the Meyer-Wallach 
global entanglement measure 29]. 

In this paper, we study the changes of entanglement for the initial four-qubit standard W 
state which will be mixed due to transmission through noisy channels. For this purpose, we 
prepare initially a pure W state then by solving analytically Eq. (Jl]) with Pauli matrices as 
the Lindblad operators we obtain dynamical evolution of system under the influence of Pauli 
channels as well as isotropic (depolarizing) channel. Also, to compute the entanglement, we 


employ a lower bound for the multi-qubit concurrence proposed by Li et al. 
that the W state with the form 


30] - We note 


\W N ) = —= (|00 ... 01) + |00 ... 10) + ... + |01... 00) + 110 ... 00)), 
v N 


( 2 ) 


is an example of multipartite entangled states in which N denotes the number of qubits in 
the state. Next, in order to compare the robustness 1 of W 4 state with maximally entangled 
four-partite states, we examine the behavior of entanglement for three four-qubit maximally 


1 The notion of “robustness” here and throughout the paper is the persistence of the lower bound of 
entanglement in noisy environments. 
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FIG. 1: A schematic quantum circuit for the transmission of a four-qubit state through a noisy 
environment. 


entangled states that are given by [31] 


l*) = ^| ( l 0000 ) + l 1111 > ) , (3) 

I&) = 1 (11111) + 11100) + |0010) + 10001)), (4) 

\h) = -j= (V2|llll) + |1000) + |0100) + |0010) + |0001)) . (5) 

The rest of paper is as follows: In the next section, we analytically solve the Lindblad 
equation for the initial four-qubit W state and study the evolution of state under some 
Markovian noises. Also, we study the behavior of entanglement of states under the Pauli 
noises. In Sec. m we investigate the time evolution of entanglement for maximally entangled 
four-qubit states. In the last section, we present our conclusions. 


II. EVOLUTION OF ENTANGLEMENT OF W 4 STATE UNDER NOISY CHAN¬ 
NELS 

In this section, we analytically solve the Lindblad equation ([Tj) for the initial W 4 state 
in contact with some Markovian noises and study the dynamics of the entanglement of the 
states. It is remarkable to point that the Lindblad equation for a four-particle state consists 
of 136 coupled differential equations which in principle their solution is a difficult task. Here, 
using the evolution of the density matrix at small times, we propose a proper ansatz for the 
sought-after density matrix which leads to at most 11 coupled equations that is much simpler 


4 















to solve. Fig. (JTJ) shows the quantum circuit that describes schematically the effect of noise 
on the state. 

For the first case, consider Pauli-X (bit-flip) channel. When a four-qubit W state trans¬ 
mits through this noise, Eq. ([I]) leads to 16 diagonal and 120 off-diagonal coupled linear 
differential equations. Thus, to simplify the problem, we first consider the time evolution of 
the density matrix W 4 state for infinitesimal time interval St as 


p(St) = p( 0) + 


' 4 

X 

i=1 


L iiX p(0)L\ 


2 i > 


St, 


in which L, T = 


and 


/ 


P{ 0) = 


Substituting p(0) in Eq. © gives 
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Now, consider the following ansatz for all times which is consistent with p^ (St): 
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( 6 ) 


(7) 


( 8 ) 


(9) 
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Inserting this solution in the Lindblad equation results is a set of 11 coupled differential 
equations 


a(t) = 4 k(b(t) — a(f)^, 
b(t ) = k (a(t) — 4 bit) + 3c(f) j, 
c(f) = 2k (d(t) + m(t) — 2 c(f) j, 
d(t ) = 2k[c{t ) + p(f) — 2 d(f)j, 
e(t) = 2k(b{t) + /(f) - 2e(f)), 

< /(f) = fc(3e(f) + /i(f) -4/(f)), (10) 

Mf) = 4A; (/(f) - h(f)), 
p(f) = A;^d(f) + 2m(f) - 4p(f) + g(f)j, 
m(f) = k[c(t) — 4 m(f) + n(f) + 2 p(f)j, 
h(f) = 2A; ^m(f) + g(f) — 2 n(f) j, 
k ?(*) = 2k(n(t) +p(t ) -2g(f)), 


subject to the initial conditions &(0) = c(0) = 1/4 and a(0) = d(0) = e(0) = /(0) = /r(0) = 
p(0) = m(0) = n(0) = g(0) = 0. The solutions read 


— 2 Kt 


2 e 


—6hx,t 


8 Kt 


), 
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—Gut 




2 e ~ AKt + e~ 8Kt ). 


( 11 ) 


. P(t) = m (1 

Now, using Pyy 4 (t) we examine the evolution of the entanglement for this state. As a quan¬ 
tifier of the entanglement, we use the lower bound to multipartite concurrence which is 


introduced by Li ef al. 


30] 


C N (p) > r N (p) 



N K 


EEk> ! . 


( 12 ) 


that it is given in terms of N bipartite concurrences C n correspond to the possible bipartite 
cuts of the multi-qubit system in which one of qubits is separated from the other qubits. 
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The bipartite concurrence C n is defined by a sum of K = 2 N ~ 2 {2 N ~ 1 — 1) terms which are 
expressed as = max{0, A/, — A*: — A| — A*;}, where A™, m = 1,..., 4, are the square roots 
of four nonvanishing eigenvalues in decreasing order of the matrix pp^. Here, p\ = S^p*S% 
and S% are defined by the tensor product of the generators of S0(2) and S0(2 W_1 ). For this 
case, the lower bound can be obtained numerically which is plotted in Fig. [2] (blue line). 

For the next case, consider the Pauli-Y (bit-phase-flip) noise. Calculations in this case 
show that results are similar to the previous case and therefore both noises have the same 
effect on the four-qubit W state, namely p y Wi (t) = p^(t). 

Now, consider the case that W 4 state is coupled to the Pauli-Z (phase-flip) channel. For 
this case, the infinitesimal time evolution of the density matrix reads 
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Inserting p^ 4 to the Lindblad equation (JT]) gives rise in 

a(t) = 0, 
b(t) = —4 kb(t), 

subject to the initial conditions a(0) = 6(0) = 1/4. Then, we obtain 


(13) 


(14) 


(15) 


a(t) = f, 

b(t) = 


The lower bound in this case is given by 


(16) 


n(p z w) = 


(17) 
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FIG. 2: The lower bound for the four-qubit concurrence for an initial W state transmitted through 


Pauli-X and Pauli-Y (blue line), Pauli-Z (red line), and isotropic (black line) channels as functions 


of tit. 


This result is depicted in Fig. [2] as a red line. 


To this end, we study the interaction of an initial four-qubit W state with the depolarizing 


channel. Taking into account the effects of all Pauli matrices as noisy operators in the 


Lindblad equation, the infinitesimal time evolution of the density matrix becomes 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

2 K&t 2 KSt 

0 

0 

2 KSt 2 K5t 

0 4 KSt 0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 




0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 / 



\ 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0/ 


Therefore, we take the following ansatz: 


/ 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

°\ 


0 

a 

b 

0 

b 

0 

0 

0 

b 

0 

0 

0 

0 

0 

0 

0 


0 

b 

a 

0 

b 

0 

0 

0 

b 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

e 

0 

c 

c 

0 

0 

c 

c 

0 

0 

0 

0 

0 


0 

b 

b 

0 

a 

0 

0 

0 

b 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

c 

0 

e 

c 

0 

0 

c 

0 

0 

c 

0 

0 

0 


0 

0 

0 

c 

0 

c 

e 

0 

0 

0 

c 

0 

c 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

h 

0 

0 

0 

d 

0 

d 

d 

0 


0 

b 

b 

0 

b 

0 

0 

0 

a 

0 

0 

0 

c 

0 

0 

0 


0 

0 

0 

c 

0 

c 

0 

0 

0 

e 

c 

0 

c 

0 

0 

0 


0 

0 

0 

c 

0 

0 

c 

0 

0 

c 

e 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

d 

0 

0 

0 

h 

0 

d 

d 

0 


0 

0 

0 

0 

0 

c 

c 

0 

0 

c 

c 

0 

e 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

d 

0 

0 

0 

d 

0 

h 

d 

0 


0 

0 

0 

0 

0 

0 

0 

d 

0 

0 

0 

d 

0 

d 

h 

0 

\ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

f) 
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Inserting this density matrix in Eq. ©, we find 


a(t) = 2 k(g(t) + 3 e(t) — 4a(t)j, 

b(t ) = 4 k[c(t) — 36(f) j, 

c(t ) = 2k (b(t) + d(t ) — 6 c(t) j , 

<i(t) = Akicit) — 3 d(t)\ 

< \ / (20) 
e(f) = 4/c^a(t) + 6.(6) — 2e(t)j, 

/(*) = Sk(h(t) - /(«)), 

g(t) = 8k(a(t) -g(t)), 

h(t ) = 2k (se(t) + f(t) — 46(f) ^, 

in which the initial conditions are a(0) = 6(0) = 1/4 and c(0) = d(0) = e(0) = /(0) = 
(]{{)) = 6(0) = 0. The solutions are found as 


a(t) 


6(f) = T—e 16nt 


(1 + e 4Kf ) 2 , 


c(f) = J6 e &Kt (1 — e Snt ), 
d(t) = ±e- WKt (-1 + e 4 ^) 2 , 

' e(f) = i (1 - e- 16 ^), 

/(f) = T (i _ 2e _4Kt + 2e _12Kt - e ~ 16Kt ), 

c/(f) = 2L (i + 2e~ 4Kt - 2e _12Kt - e" 16Kt ), 

k 6(f) = jq (1 - 2e~ AKt - 2e~ l2Kt + e" 16 "*). 

The corresponding lower bound for this state is depicted in Fig. [2] as a black line. 


( 21 ) 


III. EVOLUTION OF ENTANGLEMENT OF MAXIMAL ENTANGLED FOUR- 
QUBIT STATES 

In this section, we compare the time evolution of entanglement of W 4 state with those 
of maximally entangled four-qubit states, namely Eq. (l3lh For this purpose, we solve the 
Lindblad equation (JT]) for these states and then compute the lower bound to the concurrence 
for each case. 

First, let us consider the state |</>i) which is indeed the four-qubit GHZ state. For this 
case, the Lindblad equation is solved for the Pauli channels a x , a y , a z as well as the isotropic 
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channel and the lower bound dT2]l is computed for each case in Ref. 32], Here, we only 
present the results as follows: 


T x (<f> i) = T y ((j) i) = \[2 max 0, - (e 8Kt + 6e AKt 


— 8 K,t 


t 2 (0i) = V2 e 
Td((f>i) = V2 max 


0, g (9e 


— 16/ct 1 /?_ — 8ftt 


+ 6e 


( 22 ) 

(23) 

(24) 


Now, consider the state |0 2 ) which is transmitted through the Pauli-X channel. Its time 
evolution as a solution of the Lindblad equation (JTj) and after following the steps similar to 
the previous section is obtained as 


a 0 
0 6 6 


0 6 6 
a 0 


0 a 0 
0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 6 
a 0 
a 0 
.066 


0 0 
0 a 0 


0 0c 
0 0 0 
0 0 0 
0 0c 
0 0 0 
0 0c 
0 0c 
0 0 0 
6 0 0 
0 a 0 
0 a 0 
0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 c 
c c 0 
c c 0 
0 0c 
c c 0 
0 0c 
0 0c 
c c 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 c c 
c 0 0 
c 0 0 
0 c c 
c 0 0 
0 c c 
0 c c 
c 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0a 
0 6 0 0 


a 0 ' 


0 6 0 0 
0 0 a 
0 0 0 
c 0 0 
c 0 0 
0 0 0 
c 


0 0 
0 0 0 
0 0 0 
c 0 0 
0 6 0 0 
0 0 a 
0 0 a 
0 6 0 0 6 


6 \ 

6 

a 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
6 

a 0 
a 0 j 


(25) 


where 


a = i (1 + e~ 4Kt - 2e ~ 6Kt ), 
b = i (1 + e~ 4Kt + 2e ~ 6Kt ), 
^ c = hi (1 — e~ 4Kt ) ■ 

For this case, the lower bound (TT2T) becomes 


(26) 


r x {<j ) 2 ) = \p2 max <j 0, - (2e 6Kt + e 


—4 K,t 


(27) 


For the next case, the state \(j> 2 ) is exposed to the Pauli-Y channel. The solution reads 


/ 


41 _ 


a 0 0 

0 6 g 

0 g 6 

e 0 0 

0 0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 mm 
n 0 0 

n 0 0 

0mm 


e 0 0 
0 0 0 
0 0 0 
a 0 0 
0 c 0 
0 0c 
0 0 / 
0/0 
0 0 d 
0 d 0 
0 d 0 
0 0 d 
0 0 0 
n 0 0 
n 0 0 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0/0 
/ 0 d 
c 0 d 
0 c 0 
d 0 c 
0 d 0 
0 d 0 
dO f 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
d d 0 
0 0 d 
0 0 d 
d d 0 
0 0 / 
c / 0 
/ c 0 
0 0c 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 n n 0 \ 
m 0 0 m ' 
m 0 0 m 
0 n n 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
6 0 0 
0 a e 
0 e a 

g 0 0 


(28) 
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so that 


' a = JL (l + e ~ AKt - 2e~ 6Kt ), 
b = i (1 + e ~ AKt + 2e ~ 6Kt ), 
c = ± (1 « e ~ 4Kt ), 
d = i (e~ 6Kt - e ~ 2Kt ), 

< e = __L (2e~ 2Kt - e~ Ant - e~ 8Kt ), (29) 

/ = re ( e ' 4Kt - e ~ 8Kt ) > 
g = i (2eT 2nt + e _4ftt + e" 8Kt ), 

m = i (e _2Kt + 2e~ AKt + e" 6 ^), 

n=± (e~ 2nt - 2e~ AKt + e~ 6Kt ). 

lb v ' 


The lower bound (fT2l) now becomes 

r y ((j) 2 ) = V2 max jo, ^ (e" 8Kt + 4 e - 6Kt + 6e" 4Kt + 4e.- 2K ' - 7) j . (30) 

For transmission of \<p 2 ) through the Pauli-Z channel, the time evolution of the corre¬ 
sponding density matrix is 


2 


0 0 
0 


0 


0 0 0 
c 0 0 
a 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
b 0 0 
0 0 0 
0 0 0 


a 
c 

0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 b 
0 0 
0 0 
.0 b b 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 
0 6 0 
0 6 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 a 0 
0 0 0 
0 0 0 
0 c 0 


\ 


0 0 
0 6 
0 6 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 c 
0 0 
0 0 / 
0 a/ 


(31) 


where 


a = 


1 

4’ 

1 s> — 6K,t 


b={e 
C = \e~ AKt . 


(32) 


The lower bound to the concurrence for this case is similar to the lower bound obtained for 
Pauli-X channel, namely r z ((j) 2 ) = t x (0 2 )- 

Finally, if \(p 2 ) is exposed to the depolarizing channel, its density matrix takes the follow- 
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ing form 


where 



/ a 0 0 e 
/ 0 b d 0 
0 d b 0 
e 0 0 a 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 m m 0 
n 0 0 n 
n 0 0 n 
0 m m 0 


0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
c 0 0 / 0 0 
0 c / 0 0 0 
0 / c 0 0 0 
/ 0 0 c 0 0 
0 0 0 0 c 0 
0 0 0 0 0 c 
0 0 0 0 0 / 
0 0 0 0 /0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 


0 0 0 n n 0 \ 

0 0 m 0 0 m ' 

0 0 m 0 0 m 

0 0 0 n n 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 / 0 0 0 0 

/ 0 0 0 0 0 

c 0 0 0 0 0 

0 c 0 0 0 0 

00 HO d 
0 0 0 a e 0 

0 0 0 e a 0 ; 

OOdOO b ' 


’ a = ^ (1 + e~ 8Kt - 2e~ 12Kt ) , 
b — 1 + e~ 8nt + 2e ~ UKt ), 

c = £ (1 - e~ 8Kt ) , 

d = ^ ( e~ 8Kt + 2e~ 12Kt + e ~ WKt ), 


f _ 

i 

( 8nt 

a —16Kt\ 


16 

l e 

- e J 

m = 

1 

8 

/g—12 Kt 

| g—1 Gut 

n = 


g-12 Kt 

_ g—1 


In this case, Eq. oa gives 


r d (<fi 2 ) = V2 max jo, g (7e 16 “‘ 


+ 8e -12Kt - 



( 33 ) 


(34) 


(35) 


Now consider the state Ifo). If this state is transmitted through Pauli-X channel, the 


time evolution of its density matrix is obtained as 


/ 

a 

0 

0 

d 

0 

d 

d 

Q 

0 

d 

d 

9 

d 

<1 

Q 

°\ 


0 

g 

b 

0 

b 

0 

n 

m 

b 

0 

n 

m 

n 

m 

0 

c 


0 

b 

g 

0 

b 

n 

0 

m 

b 

n 

0 

m 

n 

0 

m 

c 


d 

0 

0 

P 

n 

f f 

0 

n 

/ / 

0 

0 

e 

e 

u 


0 

b 

b 

n 

g 

0 

0 

m 

b 

n 

n 

0 

0 

m 

m 

c 


d 

0 

n 

f 

0 

P 

/ 

0 

n 

f 

0 

e 

/ 

0 

e 

It 


d 

n 

0 

f 

0 

f 

p 

0 

n 

0 

/ 

e 

/ 

e 

0 

It 


g 

m 

m 

0 

m 

0 

0 

l 

0 

e 

e 

r 

e 

r 

r 

0 


0 

b 

b 

n 

b 

n 

n 

0 

g 

0 

0 

m 

0 

m 

m 

c 


d 

0 

n 

f 

n 

/ 

0 

e 

0 

P 

/ 

0 

/ 

0 

e 



d 

n 

0 

/ 

n 

0 

/ 

e 

0 

f 

P 

0 

/ 

e 

0 

u 


g 

m 

m 

0 

0 

e 

e 

r 

m 

0 

0 

l 

e 

r 

r 

0 


d 

n 

n 

0 

0 

/ / 

e 

0 

/ / 

e 

P 

0 

0 

u 


g 

m 

0 

e 

m 

0 

e 

r 

m 

0 

e 

r 

0 

l 

r 

0 


g 

0 

m 

e 

m 

e 

0 

r 

m 

e 

0 

r 

0 

r 

l 

0 

\ 

0 

c 

c 

u 

c 

u 

u 

0 

c 


it 

0 

u 

0 

0 

h) 


(36) 
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so that 


a = i (3 + 6e~ 4Kt - 8e~ &Kt - e ~ 8Kt ), 
b = ± (1 + 2e~ 2Kt + 2e~ AKt + 2 e - 6Kt + e" 8 ^) 

d = A (i + 2e_2K< _ 2e_6Kt - e_8K *) ’ 

/ = A (i - e_8Kt )> 

g = ^ (3 + Ae~ &Kt + e- 8Kt ), 
h = A (3 + 6e- 4 ^ + 8e~ 6Kt - e ~ 8Kt ), 
l = A (3 - 4e _6fct + e" 8 ^), 
m = A (1 - 2e _4Kt + e ~ 8Kt ), 
n = -^= (1 - e _2ft * - e" 4 ^ + e" 6Kt ). 
c = ^ (t + e _2ret + 3e _4K * + 3e ~ 6Kt ), 


(37) 


e = 


24v/2 

® — 2475 ^ 


(1 + e“ 2Kt - e" 4Kt - e ~ 6Kt ), 

,—2 K,t gg—4/^t 3 ^— 6 /st\ 


The lower bound for this density matrix is depicted in Fig. |3](a). 

The solution of the Lindblad equation when \^) is exposed Pauli-Y noise is given by 

( 


r,V — 


a 

0 

0 

—d 

0 

-d 

—d 

q' 

0 

-d 

—d 

q' 

—d 

q' 

q' 

0 \ 

0 

9 

b 

0 

b 

0 

n 

—m 

b 

0 

n 

—m 

n 

—m 

0 

d 

0 

b 

9 

0 

b 

n 

0 

—m 

b 

n 

0 

—m 

n 

0 

—m 

d 

—d 

0 

0 

P 

n 

f 

/ 

0 

n 

/ 

/ 

0 

0 

e' 

d 

u 

0 

b 

b 

n 

9 

0 

0 

—m 

b 

n 

n 

0 

0 

—m 

—m 

d 

—d 

0 

n 

f 

0 

P 

/ 

0 

n 

/ 

0 

d 

f 

0 

d 

—u 

—d 

n 

0 

f 

0 

f 

P 

0 

n 

0 

/ 

d 

f 

e' 

0 

—u 

q' 

—m 

— m 

0 

—m 

0 

0 

l 

0 

d 

d 

r 

d 

r 

r 

0 

0 

b 

b 

n 

b 

n 

n 

0 

9 

0 

0 

—m 

0 

—m 

—m 

d 

-d 

0 

n 

f 

n 

f 

0 

d 

0 

P 

/ 

0 

/ 

0 

d 

—u 

—d 

n 

0 

f 

n 

0 

/ 

d 

0 

f 

p 

0 

/ 

e' 

0 

—u 

q' 

—m 

— m 

0 

0 

d 

d 

r 

—m 

0 

0 

l 

d 

r 

r 

0 

—d 

n 

n 

0 

0 

f 

f 

d 

0 

/ 

/ 

d 

P 

0 

0 

—u 

q' 

—m 

0 

d 

—m 

0 

d 

r 

—m 

0 

d 

r 

0 

l 

r 

0 

q' 

0 

—m 

d 

—m 

d 

0 

r 

—m 

d 

0 

r 

0 

r 

l 

0 

0 

d 

d 

u 

d 

—u 

—u 

0 

d 

—u 

—u 

0 

—u 

0 

0 

h / 


(38) 


V 

where a, b, d, /, g, h, l, m, n, p, r and u are given in Eq. (1H7|) and d, e' and c( are found as 

- (3e~ 2Kt + 3e~ 4Kt + e~ 6Kt 1 -~ SKt ' 


c'= 1 


24p2 


e' = 


1 ( „ — 2 Kt _j_ g —4 K,t _ g —6 K,t ^ — 8 K,t 


’ 24 ^ 


1 ' = 5ib ( 3e_ “ 


_ 3g— 4Act _|_ g— 6ac£ g— 8ftt 


+ e 

e ~ 8nt ), 

)• 


(39) 


The lower bound for this density matrix is shown in Fig. |3Xb) . 

Now, let us consider the state |</> 3 ) which is transmitted through Pauli-Z channel. After 
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solving the Lindblad equation, the time evolution of the density matrix is obtained as 


P<t> 3 


0 0 
0 a 
0 6 
0 0 
0 b 
0 0 
0 0 
0 0 
0 b 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 c 


0 0 0 0 
6 0 6 0 
a 0 6 0 
0 0 0 0 
6 0 a 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
6 0 6 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
c 0 c 0 


0 0 0 0 0 

0 0 6 0 0 

0 0 6 0 0 

0 0 0 0 0 

0 0 6 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 a 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 c 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


\ 


0 0 
0 c 
0 c 
0 0 
0 c 
0 0 
0 0 
0 0 
0 c 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0/ 
0 df 


( 40 ) 


where 


a = 


6 ’ 


h — 

u 6 e , 

n _ 1 „—6Kt 

c ~ 3s/2 e 


d= §• 

Computation of the lower bound gives rise to 


(41) 


t 2 (0 3 ) = - v 7 (3e _12ftt + e~ 8Kt ) 


(42) 


To this end, consider the transmission of |c/> 3 ) through the depolarizing channel. In this 


case, the time evolution of the state is described by 


n d 
P<t> 3 


/ 


a 0 0 0 
0 d b 0 
0 6 d 0 
0 0 0 s 
0 6 6 0 
0 0 0 m 
0 0 0 m 
e 0 0 0 
0 6 6 0 
0 0 0 m 
0 0 0 m 
e 0 0 0 
0 0 0 0 
e 0 0 0 
e 0 0 0 
0 c c 0 


0 0 
6 0 
6 0 
0 m 
d 0 
0 9 
0 m 
0 0 
6 0 
0 m 
0 0 
0 0 
0 m 
0 0 
0 0 
c 0 


0 e 
0 0 
0 0 
m 0 
0 0 
in 0 
9 0 
0 / 
0 0 
0 0 
in 0 
0 n 
in 0 
0 n 
0 n 
0 0 


0 0 
6 0 
6 0 
0 m 
6 0 
0 m 
0 0 
0 0 
d 0 
0 9 
0 m 
0 0 
0 m 
0 0 
0 0 
c 0 


0 e 
0 0 
0 0 
m 0 
0 0 
0 0 
m 0 
0 n 
0 0 
m 0 
9 0 
0 / 
m 0 
0 n 
0 n 
0 0 


0 e 
0 0 
0 0 
0 0 
0 0 
m 0 
m 0 
0 n 
0 0 
m 0 
m 0 
0 n 
9 0 
0 / 
0 n 
0 0 


e 0 ■ 
0 c 
0 c 
0 0 
0 c 
0 0 
0 0 
n 0 
0 c 
0 0 
0 0 
n 0 
0 0 
n 0 
/ 0 
0 h ■ 


(43) 


so that 


a = i* (3 + 6e 


— 8 Kt 


gg—12 K,t _ g—1 


b = ± ( e~ 8Kt + 2e~ 12Kt + e ~ 16Kt ) 

r _ 1 ( p —Vint I 16 Kt\ 

C “ 1 + > ’ 

d=±(3 + 4e~ 12Kt + e ~ 16Kt ), 

{ e = (e~ 12Kt - e ~ 16Kt ), 


6\/2 

/ = 7s (3 - 4e- 12Kt 




— 1 6Kt\ 


h = -h (3 + 6e~ 8Kt + 8e~ 12Kt - e~ lM ) 


48 

m = ^ ( e - 8Kt — e" 16 ^) 

8 Kt _ 2g — 12fct _i_ p—1 


+ e“ 


( 44 ) 
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The lower bound to the concurrence is plotted in Fig. [3|d). 

In Fig. ((3]), we plotted the results for the four investigated states. It should be noted that 
although the states introduced in Eq. (ED are maximally entangled, the initial lower bound 
to the concurrence for these states gives rise to different values, i.e., r(</>i) = r(</> 2 ) = V2 
and t(</> 3 ) = 4/3. So, we normalized them so that r(</> 1 ) = r(</> 2 ) = t(</> 3 ) — 1 at nt — 0. As 
it is shown in Fig. (ED, the entanglement present in W 4 state is more robust than the other 
investigated states against decoherence under bit-flip, bit-phase-flip, and phase-flip noises. 
However, the four-qubit GHZ state is more robust in the presence of the depolarizing noise. 

(a) (b) 




(c) (d) 



FIG. 3: The lower bound to the concurrence for W 4 (blue line), four-qubit GHZ (black line), 0 2 
(red line) and </> 3 (brown line) under: (a) Pauli-X (b) Pauli-Y (c) Pauli-Z and (d) isotropic noisy 
channels versus nt. 
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IV. CONCLUSIONS 


In this paper, we have investigated the dynamics of an initial four-qubit W state in in¬ 
teraction with its surrounding environment. We exactly solved the Lindblad equation in 
which the Lindblad operators are proportional to the Pauli matrices and obtained the den¬ 
sity matrices corresponding to each noisy channel. We also examined the time evolution of 
entanglement using the lower bound to the concurrence for four-qubit W state. It is found 
that the entanglement of states vanishes after some finite time for the Pauli channels a x and 
Cy as well as the isotropic channel. However, for the Pauli-Z noise, the entanglement expo¬ 
nentially decreases and vanishes asymptotically. Moreover, by exactly solving the Lindblad 
equation, we studied the effects of various noises on three four-qubit maximally entangled 
states and obtained the time evolution of the lower bound. We found that except the de¬ 
polarizing channel, the W 4 state is more robust against the decoherence with respect to <f> 1 , 
</> 2 , and (j) 3 . Also, except Pauli-Z channel, W 3 is more robust than W 4 under the noises 25]. 
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